Exclusive diffractive processes in electron-ion collisions 
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We present a new technique to calculate the cross-section for diffractive vector meson production 
and DVCS in electron-ion collisions based on the dipole model. The measurement of these processes 
can provide valuable information on non-linear QCD phenomena, such as gluon saturation, and 
is the the only known way to gain insight into the spatial distribution of gluons in nuclei. We 
present predictions of differential cross-section distribution da/dQ 2 and da/dt for J/ip and <j> meson 
production for diffractive processes of heavy nuclei and demonstrate the feasibility of extracting the 
gluon source distribution of heavy nuclei, F(b), from coherent diffraction. We briefly introduce a new 
event generator based on our method that can be used for studying exclusive diffractive processes 
at a future electron-ion collider. 



I. INTRODUCTION 

The HERA accelerator at DESY, Germany, with colli- 
sion energies of y/s = 320 GeV was the hitherto highest 
energy lepton-hadron collider. One of the great achieve- 
ments of HERA was the determination of the partonic 
structure of the proton A lepton-hadron collision 
is mediated by a virtual photon, which interacts with a 
valence- or sea-quark within the hadron at a resolution 
Q 2 . When probed at higher energies, gluons fluctuating 
into gluon- or quark-pairs can be resolved at smaller time 
scales, such that more partons share the hadron's longitu- 
dinal momentum at higher energies. At small momentum 
fractions x < 1CP 2 of the participating partons, measure- 
ments at HERA showed that the content of the proton 
is dominated by gluons, and that the gluon number den- 
sity at smaller x seems to rise uncontrollably. When ex- 
trapolating current measurements to small x-values, the 
gluonic part of the cross-section becomes larger than the 
total proton cross-section. This violation of the unitar- 
ity bound can only be avoided by introducing satura- 
tion effects that tame the explosive growth of the gluon 
density. While many saturation models describing these 
non-linear effects were developed 0, H| , there exists no di- 
rect measurement that would allow to verify these models 
and ultimately prove the existence of gluon saturation. 
Although more and more tantalizing hints of the onset 
of gluon saturation coming from proton-ion collisions at 
RHIC have become available , alternative explana- 
tions can currently not be ruled out flol - [T2l ]. The direct 
study of these non-linear saturation effects would require 
lepton-hadron collisions at energies far exceeding those at 
HERA. Electron-ion collisions offer an alternative way to 
study high gluon-density phenomena at an order of mag- 
nitude lower center-of-mass energies. At high enough en- 
ergies the small x gluons in the heavy ion have a wave 
length in the longitudinal direction that encompass the 
entire width of the nucleus. A probe will thus coherently 
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interact with the bulk of low-x gluons. For a heavy ion, 
the thickness is approximately constant away from the 
edges and is proportional to A 1 / 3 , where A is its atomic 
number. This a ppr oximate dependence is supported by 
detailed studies jl3l ITij . Therefore, probing a heavy ion 
with A w 200 is similar to probing a proton at 6 times 
higher energy, making the nucleus an efficient amplifier 
of the physics of high gluon densities. 

There are two proposed future collider projects that 
aim to directly measure the saturated gluon regime 
for the first time: the Large Hadron-electron Collider 
(LHeC) at CERN using the existing LHC machine [HI 
and the Electron-Ion Collider (EIC) in the USA [l7|, us- 
ing either the existing RHIC accelerator complex at BNL 
(eRHIC), or the existing electron beams at JLab (MEIC). 

At HERA, an unexpected discovery was that approx- 
imately 10% of the ep cross-section is from diffractive 
final states [lj| and that this fraction is fairly indepen- 
dent of W and Q 2 . What characterizes these events ex- 
perimentally is the presence of a rapidity gap, a region 
in the angular coverage which exhibits no hadronic ac- 
tivity. Diffractive interactions result when the electron 
probe in DIS interacts with a color neutral vacuum exci- 
tation. This vacuum excitation, which in perturbativc 
QCD may be visualized as a colorless combination of 
two or more gluons, is often called the Pomeron. The 
hard diffractive cross-section is proportional to the gluon- 
density squared, making it the most sensitive probe of 
gluon density known. Thus, diffraction and saturation 
are closely related phenomena. 

Measurements of diffraction in an electron-ion collider 
also have substantial potential to shed light on other 
unanswered questions in heavy ion collision [l7|. Mea- 
surements over the last decade in heavy ion collision ex- 
periments at RHIC indicate the formation of a strongly 
coupled plasma of quarks and gluons (sQGP). This sQGP 
appears to behave like a "near-perfect liquid" with a ratio 
of the shear viscosity to entropy density (rj/s) approach- 
ing 1/Att [l8l - [2l| . Recent experiments at the LHC with 
substantially higher energies and thus a hotter and longer 
lived plasma phase confirm this picture [22[ ■ Despite the 
significant insight that the sQGP is a strongly correlated 
nearly perfect liquid, little is understood about how the 
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system is created. The largest uncertainty in our under- 
standing of the evolution of a heavy-ion collision comes 
from our limited knowledge of the initial condition, i.e. 
momentum and spatial distributions of gluons in the nu- 
clei. Also of importance is how the spatial distribution 
fluctuates around its mean, since it affects the behavior 
of collective effects such as flow and their higher mo- 
ments. For example, different assumptions about the nu- 
clear initial distributions give differences up to factors of 
two for the obtained r]/s value 0, [Hj|. Measurements 
of the initial gluon distribution with existing machines 
are only possible indirectly and with large uncertainties. 
The study of gluon distributions using diffractive events 
in electron-ion collisions would allow one to directly mea- 
sure the initial condition of the colliding ions, providing 
both its momentum and spatial distributions as well as 
the underlying fluctuations ("lumpiness"). In fact, ex- 
clusive diffractive eA events are the only way to study 
the initial spatial distributions and shed light on these 
fundamental questions. 

In a diffractive eA event, the electron collides with the 
ion producing one or more extra particles but leaving 
the nucleus intact. The interaction with the nucleus is 
either elastic or inelastic, and in the latter case the nu- 
cleus subsequently radiates a photon or breaks up into 
color neutral fragments. When it stays intact, the event 
is called coherent and when it breaks, the event is called 
incoherent. The spectrum of the cross-section with re- 
spect to the hadronic momentum transfer t is related to 
the transverse spatial distribution of the gluons in the 
ion through a Fourier transform. Also, according to the 
Good- Walker picture [25| , the incoherent cross-section is 
a direct measure of the lumpiness of the gluons in the 
ion. In order to access t in these events, the complete 
final state has to be measured. This is experimentally 
only possible in events such as vector-meson production 
or Deeply Virtual Compton Scattering (DVCS). 

At present, the most common approach to calculate 
diffractive cross-sections at small x is in the dipole model, 
where the exchanged virtual photon splits up into a quark 
anti-quark pair, which forms a color dipole. The dipole 
subsequently interacts with the target in the target's rest 
frame. The dipole model became an important tool for 
DIS when Golec-Biernat and Wusthoff (GBW) [H, [27| 
observed that a simple ansatz for the dipole model inte- 
grated over the impact parameter was able to simulta- 
neously describe the total inclusive and diffractive cross- 
sections. The GBW model also naturally contains sat- 
uration in the small x regime. A shortcoming of the 
GBW model is that it cannot describe the high Q 2 scal- 
ing violation in the inclusive cross-sections measured at 
HERA, something perfectly described by the collinear 
DGLAP formalism, which in turn cannot describe the 
high fraction of diffractive events. This sparked Bartels, 
Golec-Biernat and Kowalski (BGBK) to include an ex- 
plicit DGLAP gluon distribution into the dipole formal- 
ism (28|, taken at a scale directly linked to the dipole 
size. The BGBK model replicates the GBW model where 



it is applicable and also manages do describe the Q 2 de- 
pendence of the cross-sections. However, this approach 
still integrates out the impact parameter dependence of 
the interaction, without which the ^dependence of the 
cross-section is unknown. The impact parameter depen- 
dence was introduced in the dipole model by Kowalski 
and Teaney [30| and then modified to also include ex- 
clusive processes by Kowalski, Teaney and Motyka [3l| . 
This dipole model goes by the name bSat (or sometimes 
IPSat), and is the main focus of this paper. 

The bSat model has been studied in detail in the case 
of electron proton collisions at HERA. There are a few 
theoretical attempts to expand the bSat model to also de- 
scribe exclusive eA collisions (see e.g. [3(1 [HI, HH ) . With- 
out exception, these models fail to describe the disap- 
pearance of the incoherent cross-section as t — > 0. Also, 
they turn out to be poorly suited for implementation in 
a Monte Carlo event generator. 

In this paper we present the first calculations of not 
only the coherent but also the incoherent cross-sections 
in electron-ion collision without making approximations 
larger than those already inherently present in the bSat 
model, for all t. We have implemented the calculation 
described in this paper in a Monte Carlo event generator 
(Sartre). 

The paper is organized as followed: In section|H]we will 
show our derivation of the dipole model in eA, taking as 
a starting point the case of ep. In section IIIII we will 
present the resulting cross-sections, both as comparisons 
with HERA data and as predictions for EIC and RHIC. 



II. THE BSAT DIPOLE MODEL 

Earlier studies of the dipole model showed that a wide 
variety of DIS data can be described with only a few as- 
sumptions. In particular, it was demonstrated that inclu- 
sive DIS can be described together with inclusive charm 
production and exclusive diffractive vector meson photo- 
and electro-production. Especially the bSat dipole model 
is very successful in describing the exclusive production 
of J/V>, 0, p, and photon (DVCS) production at HERA. 
Here we only give a short overview of the bSat model in 
ep before we discuss its extension to eA collisions. For a 
detailed discussion on the bSat model see l3ll. 



A. A brief description of the bSat dipole model in 
diffractive ep 

The amplitude for producing an exclusive vector meson 
or a real photon diffractively in DIS can be written as: 

A££* v *(x,Q,A)=i JdrJ^J d 2 b(*^) (r,z) 

da ip) 

x27rrJ ([l - z ] r A)e~ lb - A -^-(x, r, b) (1) 
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V = J/ip, <j>, p, y 



P'.A' 



FIG. 1. A schematic picture of the dipole model and its vari- 
ables. See text for details. 



where T and L represent the transverse and longitudi- 
nal polarizations of the virtual photon, r is the size of 
the dipole, z the energy fraction of the photon taken 
by the quark, A = \f—i is the transverse part of the 
four-momentum difference of the outgoing and incom- 
ing proton, and b is the impact parameter of the dipole 
relative to the proton (see Fig. [T]). (^y^/) denote the 
wave-function overlap between the virtual photon and 
the produced vector meson. In this paper we use the 
"boosted Gaussian" wave-overlap with the parameters 
given in [3l[ . 

The dipole cross-section da^ /d 2 b(x, r, b) is defined 



da, 



Cp) 



d 2 b 



(x,r,b) = 2Af ip) {x,r,h) = 2[1 - (2) 



The first equality is the optical theorem, and we make the 
approximation of only using the real part of the S'-matrix 
for the definition of the scattering amplitude J\f ', which 
then becomes a real number between and 1. Here (p) 
denotes proton. 

In the bSat model the scattering amplitude is: 



J\fW(x,r,b) = l-e 



: r 2 a s (^)xg(x, fJ , 2 )T(b) 



(3) 



where /i 2 = 4/r 2 + /i 2 , and ^ is a cut-off scale in the 
DGLAP evolution of the gluons. The initial gluon den- 
sity xg(x,fj,Q) = A g x~ Xs (l — x) 5 - 6 . The nucleon profile 
function T(b) = l/(2irB G )cxp(-b 2 /(2B G )). All param- 
eter values are determined through fits to HERA data 
[3l| . For all results in this paper, we use Bq = 4 GcV~ 2 , 
(4 = 1-17 GeV 2 , X g = 0.02, and A g = 2.55. Also, the 
four lightest quark masses are treated as parameters in 
the model, and are taken to be: m u = m c i = m s = 0.14 
GeV, m c = 1.4 GeV. It should be noted that bSat is a 
model of multiple two-gluon exchanges at leading log, but 
some NLL effects are taken into account by the running 
of the strong coupling. 

The total diffractive j*p cross-section for this process 

is: 



dcrT P 

dt 



1 

167T 



\A(x,Q 2 ,t)\ 



(4) 



B. Extending the bSat model from ep to eA 

The explicit impact parameter dependence of the bSat 
model makes it especially well suited for the description 
of processes in eA collisions. The b dependence allows one 
to model the nucleus as a collection of nucleons according 
to a given nuclear transverse density distribution, e.g. the 
Woods-Saxon function. To this end we make two obser- 
vations. Firstly, at small x, the life-time of the dipole is so 
large that the dipole traverses the full longitudinal extent 
of the nucleus. As a consequence the nucleus can effec- 
tively be treated as a two-dimensional object in the trans- 
verse plane. Also, when the gluon's momentum fraction 
of the hadron is small, its wavelength in the light-cone 
direction x~ becomes so large, that it coherently probes 
the whole nucleus at x < A -1 / 3 / (M N R p ) ~ 10~ 2 , where 
M/v is the mass of the nucleus and R p is the proton ra- 
dius. Consequently, the information about which nucleon 
the gluon belongs to is lost, and the exact position of each 
nucleon within the nucleus is not an observable. In or- 
der to calculate the cross-section correctly the average 
over all possible states of nucleon configurations has to 
be taken: 



dcrto 



tal 



-^(\A(x,Q 2 ,t,n)\ 



(5) 



dt lb7T v ''in 
where Q denotes nucleon configurations. 

One defines two different kinds of diffractive events in 
eA: coherent and incoherent. In the Good- Walker pic- 
ture [25| the incoherent cross-section is proportional to 
the variance of the amplitude with respect to the initial 
nucleon configurations 51 of the nucleus: 



da. 



incoherent 

"dt 



-\(A(x,Q 2 ,t,n)) n 



(6) 



where the first term on the R.H.S is the total diffractive 
cross-section and the second term is the coherent part of 
the cross-section. 

When extending the bSat model from ep to eA we 
will use the independent scattering approximation to con- 
struct the scattering amplitude for nuclei: 

A 

1 - (x, r, b) = JJ (l - Af (p) (x, v, |b - b,|)) (7) 

i=i 

where b^ is the position of each nucleon in the nu- 
cleus in the transverse plane. We assume that the po- 
sitions of the nucleons are distributed according to the 
3-dimensional Woods-Saxon function projected onto the 
transverse plane. For details see Appendix lAl 

Combining equations ^ , ([3]) and ([7]) the bSat scatter- 
ing amplitude for eA becomes: 

2 d 2 b ^' ? ' b '^ ~ 



(8) 



1 — exp 



2N ( 



-r 2 a s (n 2 )xg(x,n 



2 ^E 

i=l 



T(|b-b; 
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Note that the dependence on nucleon configurations f2 in 
the amplitude is entirely contained in this dipole cross- 
section. 



1. The incoherent, coherent, and total diffractive 
cross-sections 

In order to obtain the total diffractive cross-section 
and its coherent part, the second and first moments of 
the amplitude have to be calculated respectively. For the 
first moment there is a closed expression for the average 
of the dipole cross-section [3(| : 



d 2 b 



= 2 



T A (b) p 
I 1 2~°i9 



(9) 



where is the ep dipole cross-section, eq. <j3j> , 



inte- 



grated over the impact parameter, and Ta is the profile 
of the Woods-Saxon potential in transverse space. 

For the second moment of the am plit ude, no analytical 
expression exists. Similarly as in [29(, we derive it by 
defining an average of an observable 0(f2) over nucleon 
configurations flj by: 



1 c m « 
(°>n = r~ E 





(10) 



For a large enough number of configurations C ma x the 
sum on the R.H.S. will converge to the true average. For 
the total diffractive cross-section one gets: 



d^' A 
dt 



(x,Q 2 ,t) = 



1 



1 



167rC n 



C max 

E 



A(x,Q 2 ,t,tti)\ -(11) 



For large t the variance is several orders of magnitude 
larger than the average. This means that the conver- 
gence of the sum in eq. (|10[) becomes extremely slow, as 
demonstrated in Fig. [21a), where we show the coherent 
cross-section resulting from averaging over 10, 100, 500, 
and 800 configurations. As a comparison the "analyti- 
cal average" , i.e. eq. © is also shown. As can be seen, 
not even 800 configurations are enough for convergence 
at -t > 0.15. 

The convergence of the second moment of the ampli- 
tude is shown in Fig.[2jb). We conclude that around 500 
configurations are needed to obtain a good description of 
the cross-section for —t < 0.3. 



the cross-section for small xjp ov large dipole radii in 
this case, and studies are only valid where (3 = xp/xBj 
is large. For exclusive diffraction this is equivalent to 
keeping Q 2 large. Any other way to impose a limit on the 
rise of the cross-section, e.g. through a cut-off, inevitably 
also imposes some form of saturation into the formalism. 

In the proton case, the bNonSat dipole cross-section 
is obtained by keeping the first term in the expansion of 
the exponent in the bSat dipole cross-section |30|: 



» 



-r 2 a s {p 2 )xg(x, fi 2 )T(b). 



(12) 



d 2 b N c 

In the case of a nucleus the dipole cross-section be- 
comes: 



da 



(A) 



d 2 b 



7T 

N, 



c 



r 2 a s ^ 2 )xg(x^ 2 )J2T(\b-^\) (13) 



i=i 



and the coherent part of the bNonSat cross-section can 
be obtained through the average: 

=^ 2 a s ^ 2 )xg{x^ 2 )AT A {b). (14) 

The parameters we use for the bNonSat model were ob- 
tained in [1^, by fits to HERA data. They are: B G = A 
GeV" 2 , nl = 0.8 GeV 2 , A 9 = -0.13, and A g = 3.5. The 



m d 



m s 



0.15 GeV, 



bNonSat quark masses are: 
m c = 1.4 GeV. 

Figures [3] (a) and (b) shows the wave-overlap (^>y^>) 
between the virtual photon and produced vector mesons 
as a function of dipole size r, for transverse and lon- 
gitudinal polarizations of the photon respectively. The 
wave-overlap is taken at Q 2 = 1 GeV 2 and at z = 0.7. 
In Fig. [3] (c) we show the dipole cross-section as a func- 
tion of dipole size r. In bSat the rise of the cross-section 
at large r is tamed in the model, while in bNonSat it is 
allowed to rise uncontrollably. Notice that despite the 
uncontrolled rise of the dipole cross-section, the result- 
ing cross-section stays finite because of the steep fall of 
the wave-overlap function at large r. As can be seen in 
the figure, the lighter (larger) vector mesons p and <f> are 
more sensitive to saturation effects than heavier vector 
meson such as J/tp. For J/i[> the wave-overlap falls off 
so quickly at large r that it is an unsuitable probe for 
accessing the saturated regime, even for large nuclei. 



Phenomenological corrections to the dipole cross-section 



2. A non-saturated bSat model. 

Saturation is introduced in the bSat model through the 
exponential term in the scattering amplitude (eq. ©). In 
order to study the effects of saturation on the production 
cross-section we construct a non-saturated version of the 
bSat model by linearizing the dipole cross-section. It 
should be noted that there is no taming of the rise of 



In the derivation of the dipole amplitude only the real 
part of the S'-matrix is taken into account. The imagi- 
nary part of the scattering amplitude can be included by 
multiplying the cross-section by a factor (1 + (3 2 ), where 
/3 is the ratio of the imaginary and real parts of the scat- 
tering amplitude. It is calculated using 31|: 



P = tan ( A— ) , where A 



(15) 



5 




FIG. 2. The resulting coherent (left) and total cross-section (right) for j*A — > y*J/ijjA, averaged over 10, 100, 500 and 800 
configurations. As reference, the coherent analytical average described by eq. (|9]) is also shown. 



In the derivation of the dipole amplitude, the gluons 
in the two-gluon exchange in the interaction are assumed 
to carry the same momentum fraction of the proton or 
nucleus. To take into account that they carry different 
momentum fractions, a so-called skcwedness correction is 
applied to the cross-section by multiplying it by a factor 
R g (X) defined by HJ: 

2^+ 3 T(A + 5/2) 

flg(A)= ^ r ( A + 4) (16) 

where A is defined as above. Note that this definition 
of skewedness-correction for the bSat model is slightly 
different from the one used in (3l[, but follows the de- 
scription in (33j . 

These corrections are important for describing HERA 
data, where the models are valid the corrections are typ- 
ically around 60% of the cross-section, out of which the 
skewedness correction amounts to around 15%. The cor- 
rections grow dramatically in the large x range outside 
the validity of the models, where x > 10~ 2 . 



C. Computing the eA cross-sections 

The differential ep and eA cross-sections for exclusive 
diffractive processes cannot be calculated analytically. In 
order to obtain numerical solutions we have written a 
computer program to sample and average over nuclear 
configurations. This program is also the core of a novel 
event generator, Sartre, which is briefly described in Ap- 
pendix [B] 



The total differential cross-section is: 



dQ 2 dW 2 dt 



E 

T,L 



16tt dQ 2 dW 2 



(|.4 T ,L| 2 ) f2 (17) 



where dnj L /dQ 2 dW 2 is the flux of transversely and lon- 
gitudinally polarized virtual photons, and the average 
over configurations D, is defined in eq. (|10[) . 
The coherent part of the cross-section is: 



d ^coherent 

dQ 2 dW^ 2 di 



E 

T,L 



R 2 g {l + P 2 ) dnl L 
16tt dQ 2 dW 2 



\(A 



T.LI 



'(18) 



while the incoherent part is the difference between the 
total and coherent cross-sections. 

For the the second moment of the amplitude, for each 
nucleon configuration fij, one need to calculate the inte- 
gral: 

f dz 

A t ,l(Q f,A,x P ,Qi)= I rdr—d 2 h(^) TL (Q 2 ,r,z) 

xJo([l-2]rA)e- lbA J|(.T^,r,b,a) (19) 

where the dipole cross-section is defined in eq. © for 
bSat and in eq. (fF3"]l for bNonSat. For eA, there is no an- 
gular symmetry in b which makes this integral complex. 
We average over 500 nucleon configurations, giving 1000 
such integrals for each point in phase-space. 

For the first moment of the amplitude, the integral to 
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FIG. 3. In (a) and (b) the wave-overlap between the vir- 
tual photon and produced vector mesons are shown for trans- 
verse and longitudinal polarizations respectively, as functions 
of dipole radius r. In the third panel the dipole cross-section 
is shown as a function of r, with bSat (solid) and bNonSat 
(dashed) for protons (black) and gold ions (red/grey). 



calculate is: 



xJ {[l-z]rA)J (bA) 



irrdrdzbdb (*^*) T >L {Q 2 ,r, z) 
dcr„ 



d 2 b 



(x F ,r,b)(20) 



where the average in the last term is denned in cq. (J9j) 
for bSat and in eq. ([H} for bNonSat. 

The dipole models described here are only valid for 
small values of x < 10~ 2 and not too small values of 
(3 = x/x p. If j3 becomes too small the qq dipole becomes 
unphysically large [34| . To rectify this one would need to 
include higher Fock state dipoles, such as qqg. 



III. RESULTS 

In order to verify that our numerical implementation 
reproduces measured data, we repeated the comparison 
to the latest HERA data on p, (j), J/tp, and DVCS. 
We find that both models, bSat and bNonSat, describe 
HERA data well, within the experimental uncertainties 
and within the kinematic validity of the models. This is 
not surprising since the ep part is a repetition of previous 



work ( [30|, |31| ) , although our treatment of the skewedness 
correction differs slightly. 



A. Predictions for eA collisions 

To date, there exist no experimental data on diffractive 
vector meson production in eA. However, these measure- 
ments are integral parts of the physics programs of future 
facilities such as the EIC [13 and the LHeC 0. Wc 
show results for J/ijj and </> production. We let the J/ip 
mesons decay into electron pairs, and the <fr mesons into 
kaon pairs. The pscudo rapidity and momenta of these 
decay products are restricted to \rj\ < 4 and p > 1 GeV, 
respectively. These cuts are made to limit the predictions 
to an experimentally accessible region of phase-space. We 
also limit the predictions to x < 10 -2 and Q 2 > 1 GeV 2 . 
We have simulated data corresponding to an integrated 
luminosity of 10 fb _1 , with EIC beam energies of 20 GeV 
for the electron, and 100 GcV/u for the ion beam. This 
will amount to a few months of beam operation. The 
errors shown are statistical only. 

In the left hand side panels of Figs. @] and [3] differential 
cross-sections with respect to Q 2 for J/ip and 4> produc- 
tion respectively are shown for both bSat and bNonSat 
models. The cross-sections are scaled by a factor A 4 / 3 . 
In the dilute limit (large Q 2 ) this scaling is expected to 
hold for the integral of the coherent peak, which dom- 
inates the cross-section, while deviations from it is due 
to the dense gluon regime. In the right hand side pan- 
els the ratio of ep to eAu cross-sections are shown for 
both bSat and bNonSat. As can be seen there are sig- 
nificant differences between the two models, something 
not observed at HERA. Also, the difference is larger for 
</> mesons. The reason for this is that the wave-function 
overlap between the cj> meson and virtual photon allows 
for larger dipoles than that for J ftp (sec Fig. [3]). There- 
fore, 4> production can probe further into the dense gluon 
regime and exhibits larger differences between bSat and 
bNonSat. 



Probing the spatial gluon distribution 

In Fig. |5] we show the differential cross-section with 
respect to t, da/dt, for both J/ip- and </>-meson produc- 
tion, again for both dipole models. We assume a conser- 
vative ^-resolution of 5%, which should be achievable by 
future EIC detectors. The statistical error bars shown 
correspond to an integrated luminosity of 10 fb _1 . As 
can be seen, the coherent cross-section clearly exhibits 
the typical diffractive pattern. Also depicted in Fig. [S] 
is the incoherent cross-section, which is proportional to 
the lumpiness of the nucleus. Experimentally the sum 
of the coherent and incoherent parts of the cross-section 
is measured. Through the detection of emitted neutrons 
(e.g. by zero-degree calorimeters) from the nuclear break- 
up in the incoherent case it should be experimentally fca- 
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sible to disentangle the two contributions unambiguously. 

The coherent distributions in Fig. [5] can be used to ob- 
tain information about the gluon distribution in impact- 
parameter space through a Fourier transform. In eq. (|20[) , 
the first moment of the diffractivc amplitude is a Fourier 
transform of the dipole cross-section averaged over nu- 
cleon configurations, times the wave-function overlap be- 
tween the vector meson and virtual photon. This repre- 
sents a transformation from coordinate space to momen- 
tum space A. The coherent cross-section dcr co hercnt /dt 
is proportional to the absolute square of this amplitude. 
Following {35j . we can regain the impact-parameter de- 
pendence by performing a Fourier transform on the am- 
plitude. The amplitude can be obtained by taking the 
square root of the cross-section. In order to maintain the 



oscillatory structure of the amplitude we have to switch 
its sign in every second minimum. We call this modified 
amplitude \/der co hcrcnt /dt | moc i • Its Fourier transform is: 



F(b) = 



2tt 



dAAJ (A6) 



da, 



coherent 

~dl 



(A) 



,(21) 



which is a function of impact-parameter only. In our 
models the impact-parameter dependence comes from 
the transverse density function T^(b). For bNonSat, F(b) 
is directly proportional to the input density function, 
while for bSat the relation is more complex. 

In Fig. [7] we show the resulting Fourier transforms 
of the coherent curves in Fig. [SJ using the range where 
—t < 0.36 GeV 2 . The obtained distributions have been 
normalized to unity. For testing the robustness of the 
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method, we used the statistical errors in da/dt to gener- 
ate two enveloping curves, da j dt(ti)±6 (U) , where 5 is the 
one sigma statistical error in each bin ti . The curves are 
then transformed individually, and the resulting differ- 
ence defines the uncertainty band on F(b). Surprisingly, 
the uncertainties due to the statistical error are negligi- 
ble, and are barely visible in Fig. [7] 

As a reference we show (dotted line) the original input 
distribution Ta(6), which is the Woods-Saxon function 
integrated over the longitudinal direction and normal- 



ized to unity. The bNonSat curves for </>- and J / tp-vasson 
production reproduce the shape of the input distribution 
perfectly as is expected since the bNonSat amplitude is 
directly proportional to the input distribution. For bSat, 
the shape of the J/ip curve also reproduces the input 
distribution, while the 4> curve does not. As explained 
earlier, this is not surprising, as the size of the J/ip me- 
son is much smaller than that for </>, which makes the 
latter more susceptible to differences in the dipole cross- 
section between bNonSat and bSat, as seen in Fig. El 
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We conclude that the J/ip is better suited for probing the 
transverse structure of the nucleus. However, by measur- 
ing F(b) with both J/ip- and <p mesons, one can obtain 
valuable information on how sensitive the measurement 
is to non-linear effects. Thus, both measurements are im- 
portant and complementary to each other. The results 
in Figure [7] provide a strong indication that the EIC and 
the LHeC will be able to obtain the nuclear spatial gluon 
distribution from the measured coherent t-spectrum from 
exclusive J/ip production in eA, in a model independent 
fashion. 

Strictly, the integral over A in eq. (|2Tj) should be per- 
formed up to A = oo. In Fig. [5] we demonstrate the 
effect of finite integration limits, using as an example the 
</> meson curve. We show the transformation for 4 upper 
values: |*| max = {0.025, 0.05, 0.1, 0.2} GeV 2 . The study 
shows a surprisingly fast convergence towards the input 
Woods-Saxon distribution. 



B. Ultra Peripheral Collisions 

The calculations described in this paper can also 
be applied to Ultra Peripheral Collisions (UPC) at 
hadron colliders, such as RHIC and the LHC. At very 
large impact-parameters between colliding hadrons the 
long range electromagnetic force becomes dominant over 
short-range QCD. We substitute the electron's photon 
flux dn 7 /dQ 2 dW 2 in eq. (TTT|) with that from a proton or 
an ion, as described in e.g. |36| . 



Process 


Cross-section (nb) 


P + P 


0.716 


p+Au 


0.666 • 10 a 


Au+Au 


1.22 • 10 b 



TABLE I. Cross-sections of J/ip in UPC events at RHIC. All 
cross-sections are for yj! = 200 GeV/u, 10" 6 < Q 2 < 1 GeV 2 , 
4 < W < 142 GeV, < -t < 0.3 GeV 2 . 

In Table U we list the predicted cross-sections for J/ip 
mesons produced exclusively at RHIC energy in p + p, 
p+Au, and Au+Au collisions. Each cross-section is a sum 
of the two possible photon directions in the events, such 
that symmetric beam particles are multiplied by a factor 
2, and the p+Au cross-section is the sum of the photon 
coming from the proton and from the gold-ion respec- 
tively. Especially for light mesons such as <p, these studies 
might provide new constraints for non-linear phenomena, 
such as saturation. Measurements at existing hadron col- 
liders are still limited in statistics at the time of writing 
but more detailed measurements will become available 
soon. The PHENIX experiment at RHIC measured the 
central UPC diffractive J/'i/'-production cross-section at 
V« = 200 GeV, for \r]{J/ip)\ < 0.35 corresponding to 
21 < W < 30 GeV, when the J/ip decays into an elec- 
tron pair (37j . The resulting cross-section is measured 
to be da/dy = 76 ± 33(stat.)±ll(syst.)^b. Our result is 
da/dy = 54.9 ^tb, which is well within the experimental 
uncertainty. 

IV. CONCLUSIONS 

We have presented a new method for calculating ex- 
clusive diffractive vector meson and DVCS production 
in high energy eA collisions, based on the dipole model. 
This method is the first to describe incoherent eA colli- 
sions without making approximations larger than those 
already inherently present in the dipole model, for all val- 
ues of t. In some parts of phase-space, the cross-section 
is dominated by its incoherent part, which is thus es- 
sential for making realistic predictions for future eA ex- 
periments. High energy eA collisions are expected to be 
sensitive to non-linear saturation effects. We have there- 
fore implemented our method in two dipole models: the 
bSat model and its linearization the bNonSat model. 

In Figs. H] an d O we show that in an eA collider, the 
two models are clearly distinguishable, which is not the 
case in previous ep experiments. We also show that <p- 
meson production is considerably more sensitive to non- 
linear effects than J/ip-meson production. This is due 
to the larger size of the wave-function overlap for the (p- 
meson. In Figs. [H] and [71 we show that one can probe 
the transverse spatial gluon distribution of a nucleus by 
performing a Fourier transform of the measured coher- 
ent i-spectrum. This method is very robust with respect 
to statistical uncertainties and only requires a range of 
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£ < 0.2 GeV 2 for gold. Due to its smaller wave func- 
tion, the J/i/>-meson is considerably more suitable for 
probing the spatial gluon distribution than the lighter </>- 
meson. In Fig.|5]we also show the incoherent i-spectrum, 
which is directly proportional to the lumpiness of the ini- 
tial nucleus. Our method can also be used to calculate 
UPC events in present hadron-hadron colliders. We de- 
scribe central J/tp data from the PHENIX experiment 
well within the experimental uncertainties. 



ACKNOWLEDGMENTS 

The authors would like to thank Henri Kowalski, Tuo- 
mas Lappi, and Raju Vcnugopalan for their input and 
help, and the Open Science Grid consortium for provid- 
ing resources and support. This work was supported by 
the U.S. Department of Energy under Grant No. DE- 
AC02-98CH10886. 



Appendix A: Generating a nucleon configuration 
according to the Woods-Saxon potential 

We generate the nucleus according to the Woods-Saxon 
distribution, which is assumed to describe the number 
density of nucleons per volume element, i.e.: 



d 3 iV 

~¥7 



p(r) = 



Po 



1 



(Al) 



where po is the central density, Rq is the radius of the 
nucleus and d is the skin thickness which describes how 
fast the potential falls off close to the edge of the nucleus. 
The parameters po, Rq and d have been measured for 
most nuclei in low-eneregy experiments [38[ . 

Our method for generating a nucleus is as follows: 

1. We first generate the radial distribution of all nu- 
cleons in a given nucleus specimen according to: 



AN 
dr 



and sort them in 



Airr 2 p(r) 



(A2) 



2. We then generate the angular distributions uni- 
formly in azimuthal angle, and polar angle, cos 9, 
one at a time beginning with the innermost nu- 
cleon. 

3. If the newly generated nucleon position is within 
a core distance of 0.8 fm from any other nucleon 



we regenerate <f> an d cos keeping the original r. 
If this fails repeatedly, we drop the nucleus and 
restart from 1. 

4. Finally, when all nucleons have been placed, the 
origin of the nucleus is shifted to its center of mass. 
Appendix B: Generating events with Sartre 

Sartre is a novel Monte Carlo event generator, imple- 
menting the models described in this paper. It generates 
exclusive events in diffractive vector meson and DVCS 
production for ep and eA collisions. 

The master equation of Sartre is eq. (|17[) . In the event 
generator, this cross-section is simply used as a proba- 
bility density function from which a phase-space point 
in Q 2 , W 2 , and t is drawn. Given the beam energies 
and these three kinematic variables, the final state of the 
event is fully defined except for the azimuthal angle of 
the vector meson, which is uniformly distributed. 

To determine the total cross-section in eA, the com- 
plex four-dimensional integral described in eq. (|19[) has 
to be calculated for each phase-space point 1000 times, 
which is prohibitive for efficient event generation. There- 
fore, we tabulate the first and second moments of the 
amplitudes, for both longitudinally and transversely po- 
larized photons separately. The resulting look-up tables 
are three dimensional in Q 2 , W 2 and t. There is a set of 
four look-up tables ((|-4 T | 2 ), | (At) |, (|-4l| 2 ), | (A l ) |) 
for each species of produced vector meson or DVCS, and 
for each species of nucleus. 

When an event has been generated it is decided proba- 
bilistically weather the event was coherent or incoherent 
by comparing the coherent cross-section in eq. (fT8|) with 
the total one. In the incoherent case we let the nucleus 
break up by assuming that the diffractive mass My is 
distributed according to: 



AN 



1 



AM 2 



2 K M 2 



(Bl) 



Note that My cannot be uniquely determined from kine- 
matics alone. The corresponding excitation energy of the 
nucleus is: 



E* 



(M Y -m n )-A 



(B2) 



We then use this excitation energy as input for 
Gemini+-|- (39|, a statistical model code which describes 
the nuclear de-excitation, providing the break-up prod- 
ucts from neutrons up to the heaviest fragments. 
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